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Abstract 

The nuclear shadowing and antishadowing effects are explained by a unitarized 
BFKL equation. The Q^- and 2;-variations of the nuclear parton distributions are 
detailed based on the level of the unintegrated gluon distribution. In particular, the 
asymptotical behavior of the unintegrated gluon distribution near the saturation 
limit in nuclear targets is studied. Our results in the nuclear targets are insensitive 
to the input distributions if the parameters are fixed by the data of a free proton. 
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1 Introduction 

The gluon distribution in the nuclear target is an essential ingredient in the calculation 
of high energy nuclear collisions, which relate to minijet production, dilcpton production, 
heavy quarks and their bound states. Although the gluon density in a free nucleon was 
extracted by the experiments, however there are not enough data for the gluon distribu- 
tions in nuclei. As well known the parton densities in a bound nucleon differ from in a 
free nucleon, the ratio of nuclear and deuterium structure functions is smaller or lager 
than unity at Bjorken variable x < 0.05 — 0.1 or x ~ 0.1 — 0.2. These two facts are 
called the nuclear shadowing and antishadowing effects, respectively [1]. The understand- 
ing of nuclear shadowing and antishadowing in QCD is therefore an important issue to 
predict the nuclear gluon distributions. In this aspect, the parton recombination (fusion) 
between two different nucleons in a nucleus is a natural mechanism, which transfers the 
partons from a smaller x region to a larger x region and forms the nuclear shadowing and 
antishadowing effects [2]. 

The shadowing and antishadowing phenomena are also predicted to happen in a free 
nucleon: a rapid rise of parton multiplicities inside the proton at small x leads to the 
gluon recombination, which changes the distributions of gluon- and quark-densities, but 
does not change their total momenta. In consequence, a part of the gluon momentum 
lost in the shadowing should be compensated in terms of new gluons with larger x, which 
form the antishadowing. The modification of the gluon recombination to the standard 
DGLAP [3] evolution equation was first proposed by Gribov-Levin-Ryskin and Mueller- 
Qiu (the GLR-MQ equation) in [4], and it is naturally regarded as the QCD dynamics 
of the nuclear shadowing since the type of gluon fusion in a nucleon is still there in a 
nucleus but roughly differs by an A^^^ scale [5]. However, the GLR-MQ equation does not 
predict the nuclear antishadowing effect since the momentum conservation in the gluon 



recombination is violated in this equation. For this sake, a modified DGLAP equation 
was proposed to replace the GLR-MQ equation in [6,7], where the corrections of the 
gluon fusion to the DGLAP equation lead to the shadowing and antishadowing effects. 
Unfortunately, the integral solutions either in the GLR-MQ equation or in other modified 
DGLAP equations need the initial distributions on a boundary line {x, Ql) at fixed Qq, 
and they contain the unknown input nuclear shadowing and antishadowing effects. Close, 
Qiu and Roberts [8] constructed a parton fusion model by QCD arguments and try to 
evolve the input parton distributions in the nuclear target, however, the results of this 
parton fusion model in the small x region are sensitive to the input parton distributions 
[9], which leads to a large uncertainty in the predictions. On the other hand, recently 
nuclear unintegrated gluon distribution becomes a useful and intuitive phenomenological 
language for applications to many high-energy nuclear collisions. However, the DGLAP 
equation and its modified forms are based on the coUinear factorization scheme and they 
don't predict the evolution of the unintegrated parton distributions. 

Instead of the above-mentioned modified DGLAP equations, an alternative QCD re- 
search for the nuclear shadowing and antishadowing is to modify the BFKL equation [10], 
which is directly written for the unintegrated gluon distribution. The evolutions of the 
BFKL-type equations are along the small x-direction: the input gluons distribute on the 
boundary line {xq,^) at fixed xq, where the gluon fusion begins, and all parton distribu- 
tions both in a free nucleon and in a bound nucleon x < Xq are the evolution results of 
these equations. There are several nonlinear evolution equations considering the correc- 
tions of the gluon fusion to the BFKL evolution. One of the most widely studied models 
is the Bahtsky-Kovchegov (BK) equation [11]. The BK equation is originally derived for 
the scattering amplitude in the transverse coordinator space. The nonlinear terms in the 
BK equation are formed by the dipole splitting and the screening effect origins from the 



double scattering of the probe on the target. A remarkable solution of the BK equation is 
the so-called saturation, where the amplitude is a completely flat spectrum. However, the 
amplitude in the BK equation is re-normalized to identify the coefficient of the nonlinear 
term with that of the linear term. Therefore, the relation of the unintegrated gluon distri- 
bution with this scattering amplitude is unclear and it is model-dependent [12]. Besides, 
it is similar to the GLR-MQ equation, the BK equation is irrelevant to the antishadowing 
effect. Recently, a new modified BFKL equation incorporating the shadowing and anti- 
shadowing corrections of the gluon recombination to the unintegrated gluon distribution 
was proposed by Ruan, Shen, Yang and Zhu (the RSYZ equation) [13]. 

The purpose of this work is to explain the observed nuclear shadowing and antishadow- 
ing effects for the parton distributions using the RSYZ equation. The Q^- and x- variations 
of the nuclear parton distributions in the RSYZ equation are detailed in this work. Partic- 
ularly, we predict the unintegrated nuclear gluon distributions using the RSYZ equation. 
The results show a logarithmic increasing spectral height, which is not identical to the 
prediction of the BK equation in the transverse coordinator space but is similar to a mean 
field result. 

The paper is organized as follows. In Section 2 we give the RSYZ equation and its 
modifications in the nuclear target. The predictions of the RSYZ equation to the nuclear 
shadowing and antishadowing effects in the quark- and gluon-distributions as well as in 
the unintegrated gluon distributions are presented in Section 3. 



2 The RSYZ evolution equation 

The unintegrated gluon distribution /jv in a proton obeys the RSYZ evolution equation 
at small x [13], where the contributions of the gluon recombination to the BFKL dynamics 
are considered 

dx 
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where i? = 1 /m (or < 1 /m) if the gluons are uniformly distributed in a nucleon (or the 
gluons are located in the hot-spots); xq is the starting point of the gluon fusions. Note 
that the shadowing and antishadowing coexist in the region x < xq, while there is only 
the antishadowing in xq < x < 2xo [14]. 

We should point out that the nonlinear terms in Eq. (1) are really from the evolution 
kernel of the MD-DGLAP equation [6], where the double logarithmic approximation is 
taken, i.e., both the x and transverse momenta are strongly ordered. Obviously, this 
approximation satisfies such x-region, where the values of x are not extra small. On the 
other hand, the linear terms of Eq. (1) are the BFKL-kernel, which works in the small 
X range. The mix of two approximations in Eq. (1) is feasible for the discussions of 
the nuclear shadowing and antishadowing effects, since they occupy the transition range 



from a middle x to small x. The RSYZ is directly derived for the unintegrated gluon 
distribution, which relates to the (integrated) gluon distribution using 



G^ix,Q') ^ xgr,ix,Q') = I ^fNix,k'). (2) 

Since the density of gluons increases rapidly with decreasing x, the sea quark distri- 
butions are increasingly dominated by the gluon distribution, via the DGLAP splitting 
G — > qq. Thus, the deep inelastic structure function at small x reads [15] 

F2n{x,Q') 

where Pqoixi) is the DGLAP splitting function. 

Now we discuss the RSYZ equation in the nuclear target. The gluons with smaller x 
exceed the longitudinal size of nucleon in a nucleus. We assume that the gluons inside the 
nucleus are completely overlapping and fusion along the longitudinal momentum direction 
at the evolution starting point Xq. Thus, the strength of the nonlinear recombination terms 
in Eq. (1) should be scaled by A^^^ in a nucleus. On the other hand, although the softer 
gluons of different nucleons with extra small k^ may be correlated on the transverse plane 
because the integrations on k^ can go down to a very small value in Eq. (1), we neglect 
these corrections due to fN{x,k^ ^ 0) ^ 0. In this simple model the RSYZ equation in 
the nucleus becomes 

dx 
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The distributions Ga{x, Q^) and F2a{x, Q^) in the nucleus are computed by using the 
equations corresponding to Eqs. (2) and (3). 



3 Numerical analysis and summary 

We use a parameter form of a BFKL-like solution as the input distribution [13] at 
2xq = 0.3 

U2x,£) = f42xo,e) = P^e^j—exp ) , (5) 

where Xbfkl = 12as/(7rln2) and (3 and A' are two parameters. 

The computations of the RSYZ equation need pre-know the value of fN{A){xi/2,^) 
at the step x = Xi. For this sake, we proposed the following program in [13] 
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where fN{A),Shadowing{xi/2,k^) (or fN(A),BFKL{xi/2,k^)) indicates that the evolution from 
Xi to Xi/2 is controlled by Eq. (1) but without the antishadowing contributions (or 
is controlled by the BFKL equation). The parameter rj implies the different velocities 
approaching the BK dynamics. We temporarily take rj — oo and we will indicate it is 
appropriate. 

At first, we use the well known F2n{x, (5^)-data [16] of a free proton to determine the 
parameters in the computations. Then we predict the distributions in nuclei. In this work 
we fix the coupling constant to be = 0.3. The dashed curve in Fig. 1 is our fitting 
result using P = 7.22, A' = 0.002 and R = 2.6 GeV'^. Note that the contributions from 
the valence quarks to F2 at x > 0.1 are necessary and they can be parameterized by the 
difference between the dashed curve and the experimental solid curve in Fig. 1. 

Figure 2 shows the predictions of the RSYZ equation for the Ca/C, Ca/Li, Ca/D and 
Cu/D compared with the EMC and NMC results [17,18]. The agreement is acceptable. 



Figure 3 indicates that the enhancement in the Sn gluon distribution with respect to that 
in C observed by NMC [19] is consistent with our predictions. 

No significant Q^-dependence on the ratios of the structure functions at small x has 
been concluded in the present experimental precision. However, it does not prevent us 
exposing the possible Q^-variations of the nuclear shadowing, which may be hidden in the 
larger experimental errors. The ratios Gca{x,Q^)/GD{x,Q^) for gluon distributions at 

— 2 and 10 GeV'^ using Eqs. (1) and (2) are given in Fig.4a. The Q^-variations of the 
gluon ratios are predicted in the region 10""^ < x < 10~^ in our model. The logarithmic 
slope b in Ga/ Ga' = a + & In is positive. However, the corresponding slope in the ratio 
of the structure functions F2Cal F2D is negative (see Fig.4b). For example, the predicted 
(5^-slope for calcium at x ~ 4 x 10~^, h ~ —0.046, and at x ~ 10~^, h ~ —0.03, the results 
are compatible with the measured data in [20] . A more significant Q^-dependence of the 
structure function ratios can be found in the heavy nucleus F2Xe/ F2d- A flatter ratio at 
X < 10~^ in the E665 data [21] was presented. However, it can't be understood as the 
saturation behavior, where the parton fusion balances the parton splitting. The reason 
is that the parton densities in the proton still increase toward the small x-direction at 
X < 0.01 (see Fig.l). Alternatively, we consider that this behavior is a consequence of 
the Q^-variations of the structure function ratios since the experimental point with the 
smaller x corresponds to the smaller value of in Fig. 5. 

An important parameter in the computing RSYZ equation is rj in Eq. (5). The 
observed nuclear shadowing and antishadowing provide an example to determine the 
value of T). The results incline to a minimum antishadowing, i. e., 77 = 00. 

Comparing with the modified DGLAP equation, the RSYZ equation directly predicts 
the nuclear unintegrared gluon distributions, which are important information for the 
researches of high energy nuclear collisions. We compute the unintegrated gluon distribu- 



tions in different nuclear targets using tlie RSYZ equation. In Fig. 6 we plot the spectra 
of the gluon density on the transverse momentum squared dG/d^ = dG/dQ'^\Q2^^2 — 
f{x,k^)/^ ~ ^. The BK equation on the transverse coordinator space predict the sat- 
uration solutions, where the amplitude is a completely fiat spectrum at ^ < Ql{x), and 
the spectral height is irrelevant to the values of x (i.e., the geometric scaling); Ql{x) 
is the saturation scale. The saturation limit and the geometric scahng are expected to 
appear in the larger nucleus, where the shadowing effects are enhanced. However, the 
logarithmic increasing spectral height in our results is not identical to this prediction of 
the BK equation but is similar to the mean field result: dG/dk^ oc ]n{Ql{x)/k^) [22]. 

Finally, we discuss the dependence of our solutions with the input distributions. For 
this sake, we change the parameters in Eq. (24) to: P = 0.15, A' = 9.64, R = 2.4 GeV~^ 
and strongly distort the input form. Then we repeat our calculating programs. We find 
that the results are insensitive to the input distributions. 

In summary, the nuclear shadowing and antishadowing effects are explained by a uni- 
tarized BFKL equation. The Q^- and x-variations of the nuclear parton distributions 
are detailed based on the level of the unintegrated gluon distribution. In particular, the 
asymptotical behavior of the unintegrated gluon distribution in various nuclear targets is 
studied. We find that the geometric scaling in the expected saturation range is violated. 
Our results in the nuclear targets are insensitive to the input distributions if the param- 
eters are fixed by the data of a free proton. We believe that these results are useful in 
studying the ultrarelastivistic heavy ion collisions. 
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Figure Captions 

Fig. 1 The fit of the computed F2p{x, — 10 GeV'^) in proton by the RSYZ equation 
using the input Eq. (5) (dashed curve). The contributions of the valence quarks are 
parameterized by the differences between solid and dashed curves. The data are taken 
from [16]. 

Fig. 2 Predictions of the RSYZ equation compared with the ratios of the structure 
functions for various nuclei. The data are taken from [17,18]. All curves are for Q'^ — 

10 GeV^ 

Fig. 3 Predictions for the ratio of the gluon distributions in Sn/C. The curve is the 
result of the RSYZ equation for = 10 GeV^ and the data are taken from [19]. 

Fig. 4 (a) (5^-dependence of the ratios for the gluon distributions Ca/D in the RSYZ 
equation, (b) Similar to (a) but for the the structure functions. 

Fig. 5 Q^-dependence of the ratios for the structure functions Xe/D in the RSYZ 
equation. The data are taken from [21]. 

Fig. 6 Predicted spectrum dC/dh^ = f{x,^)/^ on the fc^-space in various nuclear 
targets in the RSYZ equation. Dotted curves x — 10~^, dashed curves x — 10~^ and solid 
curves x = 10~^. 



